Continuous-time random walk with a superheavy-tailed distribution of waiting times 
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We study the long-time behavior of the probabihty density associated with the decoupled 
continuous-time random walk which is characterized by a superheavy-tailed distribution of wait- 
ing times. It is shown that if the random walk is unbiased (biased) and the jump distribution has a 
finite second moment then the properly scaled probability density converges in the long-time limit 
to a symmetric two-sided (an asymmetric one-sided) exponential density. The convergence occurs in 
such a way that all the moments of the probability density grow slower than any power of time. As 
a consequence, the reference random walk can be viewed as a generic model of superslow diffusion. 
A few examples of superheavy-tailed distributions of waiting times that give rise to qualitatively 
different laws of superslow diffusion are considered. 
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I. INTRODUCTION 

Almost a half century ago, Montroll and Weiss [l| in- 
troduced a special class of cumulative processes, the so- 
called continuous-time random walks (CTRWs). Due to 
its simplicity and flexibility, the CTRW approach has 
become an important tool for the analysis of various 
stochastic systems. It is especially eflScient for studying 
the systems whose average behavior differs from the clas- 
sical one characterized by the normal, linear dependence 
of the mean-square displacement (variance) on time. In 
particular, the CTRW approach is extensively used to 
determine the laws of anomalous diffusion occurring in 
these systems (see, e.g., Refs. 

One of the most important characteristics of the 
CTRW is the long-time behavior of the probability den- 
sity P{x, t) associated with the position x{t) of a walking 
particle. It is this density that describes in full detail the 
diffusion properties of particles, including the asymptotic 
behavior of the variance (i.e., the diffusion law). More- 
over, the probability density P{x, t) ait ^ oo plays a sig- 
nificant role in establishing the connection between the 
CTRW and the fractional time diffusion equation 5] . 

Within the CTRW framework, the probability density 
P{x,t) of the particle position is completely character- 
ized by the joint probability density of the waiting times 
Tn {n = 1,2,.. .), i.e., random times between the particle 
jumps, and the jump magnitudes a;„. For the decoupled 
CTRW, when the sets {Tn} and {xn} are statistically 
independent of each other, P{x, t) depends only on the 
probability densities p{t) and w{x) of r„ and Xn, respec- 
tively. In this case the long-time behavior of P{x, t) is 
determined by the main characteristics of p(r) and w{x) 
P-Q. For example, if the second moment of w{x) exists 
and p{t) at r — >■ cxd is given by the asymptotic formula 
p{t) ^ h/r^^" with h being a positive parameter and 
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< a < 1, then the properly scaled probability den- 
sity P{x, i) at i — >■ oo converges either to a symmetric 
two-sided (for an unbiased walk) or to an asymmetric 
one-sided (for a biased walk) stable density [6]. 

A more general class of waiting-time probability den- 
sities is described by the asymptotic relation p{t) ^ 
h{T)/T^~^°' (r — oo), where a positive function /i(r) is 
assumed to be slowly varying at infinity. Due to this 
property the function ft.(r) does not strongly influence 
the long-time behavior of x{t) if a > Q. However, at 
a — the situation changes drastically. From a mathe- 
matical point of view, this is a consequence of the fact 
that the waiting-time densities characterized by the con- 
dition a = belong to a class of superheavy-tailed densi- 
ties [13, [HI whose all fractional moments are infinite. It 
is expected, therefore, that in this case the long-time be- 
havior of the CTRW is mainly controlled by the asymp- 
totic form of slowly varying function /i(t) at r — )> oo. 
The laws of both biased and unbiased superslow diffu- 
sion fl^, i.e., diffusion of particles whose mean-square 
displacement grows slower than any power of time, cor- 
roborate this statement. 

Some important features of the long-time behavior of 
the probability density P(x, t) have already been estab- 
lished by Havlin and Weiss [Ij] for a very special class of 
slowly varying functions with /i(t) ^ a/lv}'^'^ t (a > 0, 
V > Q) as r — > oo. They have shown, in particular, that 
if the CTRW is decoupled and unbiased and the second 
moment of the jump probability density w{x) is finite 
then P(x, t) at i — > oo has symmetric exponential tails. 
In contrast, in the present paper we provide a detailed 
analysis of the long-time behavior of P{x, t) for a whole 
class of slowly varying functions /i(t) that support the 
normalization condition of p(t) and consider both biased 
and unbiased versions of the decoupled CTRW. 

The paper is organized as follows. In Sec. |TT1 we re- 
produce the basic equations of the decoupled CTRW and 
specify a class of superheavy-tailed densities of waiting 
times. The main results of the paper are obtained in 
Sec. mil Here, we determine the long-time behavior of 
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the probabihty density of the particle position and estab- 
lish its limiting form for a properly scaled spatial variable. 
Moreover, we use the two-sided exponential jump density 
to verify our results and consider a few illustrative exam- 
ples of superheavy-tailed densities of waiting times. In 
SeclIV) we calculate the moments of the particle position 
at long times and show that they grow slower than any 
power of time. In the same section, we derive the most 
general form of the laws of superslow diffusion and calcu- 
late the laws that correspond to the illustrative examples. 
Our findings are summarized in Sec. [V] 



II. SPECIFICATION OF THE MODEL 

The CTRW is defined as the cumulative continuous- 
time jump process 0-01 



where S{x) is the Dirac S function and 



N(t) 

X(t) = Y.'^r. 



(2.1) 



where X{t) (^(0) — 0) can be considered as the walker 
(particle) position, Xn € (—00,00) is the random magni- 
tude of the nth jump of a particle, and N{t) = 0, 1, 2, . . . 
is the random number of jumps up to time t (if N{t) = 
then X{t) — 0). It is assumed that both the magnitudes 
Xn and the waiting times t„ (i.e., times between succes- 
sive jumps) are independent random variables distributed 
with the probability densities w(x) andp(r), respectively. 
In the decoupled case, when the sets {xn} and {r„} of 
the variables Xn and t„ are independent of each other, 
the distribution of particles is completely determined by 
these densities. Specifically, in the Fourier-Laplace space 
the probability density P{x,t) of the particle position 
X{t) is given by the MontroU- Weiss equation [l[ 



Pk 



1 -Ps 



s(l -psWk) 



(2.2) 



Here, the indexes k and s denote the Fourier and 
Laplace transforms, Uk = F{u[x)} = J^^dxe^''^u{x) 
and V, = C{v{t)} = dte-'^vit) (Res > 0), of 
the corresponding functions. Applying to Eq. (12. 2p 
the inverse Fourier and Laplace transforms defined as 
u{x) = J^-^{uk} = (l/27r) ^^^dke-'^'^Uk and v{t) = 

C~^{vs} = (l/27ri) J^^i^ dse^*'Vs (the real parameter c 
is assumed to be larger than the real parts of all singu- 
larities of Ws), respectively, one obtains 



p{x,t)^c-'l^-^:F-' 



1 



1 I r (2.3) 

1 - PsWk J J 

It can also be shown [T3 - [l6j that the probability den- 
sity P{x, t) satisfies the integral equation 



P{x,t) ^ V{t)S{x)+ / dr / dx'P{x',T) 
Jo J -00 

X p{t — t)w{x — x'), 



V{t) = 1-1 drpir) = / drpir) 
Jo Jt 



(2.5) 



{V{0) = 1, V{oo) = 0) is the complementary cumulative 
distribution function of waiting times, that is also known 
as the survival or exceedance probability. The first term 
in the right-hand side of this equation corresponds to the 
situation, which is realized with probability V{t), when 
there are no jumps up to time t. While the exact solu- 
tion of Eq. (|2.4I) was obtained for a special case of jump 
and waiting-time densities [13], the long-time behavior 
of P(x, t) was analyzed for a much wider class of these 
densities. This has been done by using the representa- 
tion (|2.3p [1, 0| and by applying the limit theorems of 
probability theory [l^ . The first approach is based on 
the Tauberian theorem for the Laplace transforms (see, 
e.g., Ref. [T^), which states that if the function v{t) is 
ultimately monotone and 



-4- 



(0 < 7 < 00) as s — > then 
v{t) 



r(7) 



L{t) 



(2.6) 



(2.7) 



'AS t ^ 00. Here, r(7) is the gamma function and L{t) 
is a slowly varying function at infinity, i.e., L{fj,t) ~ L{t) 
as t — 00 for all > 0. We note that, in contrast to 
the inverse Laplace transform where the parameter s is 
complex with Res > 0, in Eq. ()2.6|) this parameter is 
assumed to be real and positive. 

In this paper, we use the Laplace transform 



Ps{x) 



1 -Ps 



Six) 



(1 -Ps)Ps 



Wk 



1 - PsWk 



(2.8) 

which follows from Eq. ()2.3p , and the Tauberian theorem 
to study the long-time behavior of the probability den- 
sity P{x, t) in the case of superheavy-tailed distributions 
of waiting times. More precisely, we consider a class of 
probability distributions whose densities have the follow- 
ing asymptotic behavior: 

h{r) 



p{t) 



(2.9) 



(2.4) 



(r — 00), where a positive function /i(t) varies slowly 
at infinity. Since for each slowly varying function /i(t) 
the condition r''/i(r) 00 (p > 0) holds as r — 00 
po! ]. all fractional moments of p{t) are infinite [i.e., 
Tp = dTTPpir) =00]. It should be noted that the 
slowly varying function /i(t) in Eq. (12. 9p is not arbi- 
trary: it must be compatible with the normalization con- 
dition To = /q drp^r) = 1, which in turn implies that 
/i(r) = 0(1/ In r) as r ^ 00. Finally, we assume here 
that the first two moments of the jump density wlx), 
li = J_^dxxw{x) and I2 = J_^dxx'^w{x), exist. As 
will be shown below, in this case only these characteris- 
tics of z«(x) influence the long-time behavior of P{x,t). 
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III. LONG-TIME BEHAVIOR OF P{x,t) 

According to the Tauberian theorem, the long-time be- 
havior of the probabihty density P{x, t) of the particle 
position is determined by the asymptotic behavior of the 
Laplace transform Ps{x), see Eq. (|2.8p . at s — > 0. In or- 
der to find Ps(x) for small s, let us first determine the 
two terms of the asymptotic expansion of the Laplace 
transform = dTe~^'^p{T). With this purpose, we 
use the exceedance probability (j2.5|) to represent ps as 



\s 



(3.1) 



An important feature of the exceedance probability, 
which follows from the asymptotic formula (|2.9p , is that it 
is a slowly varying function at infinity, i.e., V(iit) ^ V{t) 
as i — >■ oo [T2l|. Therefore, at s — >■ the function V{q/ s) 
in Eq. p.ip can be replaced by F(l/s), yielding 



1 - V 



(3.2) 



Since ps and Wk tend to 1 as s and k tend to 0, at small 
s the main contribution to the inverse Fourier transform 
in Eq. (|2.8p comes from the close vicinity of the point 

= 0. Due to this fact, Wk can be approximated by a 
few terms of the expansion oi Wk as fc 0. Taking into 
account that the first two moments of 'w{x), li and ^2, 
are assumed to exist, we obtain 



1 2 

Wk ^ 1 + ihk - - hk 



(3.3) 



(fc 0). Then, using Eqs. (|3?2|) and ([3^ . we can find 
the asymptotic behavior of {wk / {I ~ PsWk)} as s — > 
by replacing Wk/il—PsWk) by [hk'^ /2 + V{l/s)]~^ in the 
unbiased case (when li = 0) and by [—ilik + V{l/s)]~^ 
in the biased case (when li ^ 0). 
For the former case one has 1211 



Wk 



1 - PsWk 



1 



hk'^ 12 + Vigils) 



-\x\^2V{\ls)ll^ 



(3.4) 



According to this result, the first term in the right-hand 
side of Eq. ()2.8|) can be neglected in the limit s — 0. 
Indeed, the integral of this term over some interval con- 
taining the point x = is equal to (1 — Ps)/s ~ ^ (l/s)/s, 
while Eq. p.4|) shows that the integral of the second term 
is proportional to 1/s. Therefore, since ViXjs) — > as 
s — >■ 0, the main term of the asymptotic expansion of 
Ps{x) takes the form 



Ps(x) 



1 / nV^) |xk/2V(l/s)/i. 



2/2 



(3.5) 



(s — >■ 0). Because the exceedance probability V(€) varies 
slowly at infinity, it is not difficult to show that so also 
does the function v/V'(^)/2^2 e^l^lV^^W/'^. Hence, the 
Tauberian theorem is applicable to this case and, in ac- 
cordance with Eq. p.Sp . leads to the following behavior 
of P{x,t) at long times: 
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where, in accordance with Eq. 



y(t) 



dT-hir). 



(3.6) 



(3.7) 



It should be stressed that in spite of the fact that 
the asymptotic formula p.6p is derived considering a 
small vicinity of the point k — 0, it is valid for all x. 
The reason is that it is this vicinity which is respon- 
sible for the asymptotic behavior of Ps{x) as s — ?> 0. 
We note in this context that the asymptotic formula 
for P{x, t) obtained in [l^ for a special case of the ex- 
ceedance probability is applicable for all a;, not only for 
|a;| — > 00 as it was assumed in Ref. [isj . Finally, intro- 
ducing the variable y = xy^2V{t)/l2, we make sure that 
the limiting distribution of the scaled particle position 
Y{t) — X it) 2V {t) / 12 is described by the symmetric 
two-sided exponential density 



'P{y) = lim 



2V{t) 




2V{t) 



(3.8) 



For the biased CTRW, when h ^ 0, we obtain j2l| 



Wk 



1 - PsWk 



1 



-i/i/c + y(l/s) 



= — e-"^(i/")/'i0(/ia;), (3.9) 

where d{x) = if x < and Q{x) 1 if x > 0. It is 
important to emphasize that the condition T^^{wkl(X ~ 
VsWk)\ ~ 0, which occurs at l\x < 0, is a direct con- 
sequence of the used approximation [see the first line of 
Eq. dSl])]. In fact, the term \-ilik -h F(l/s)]-i is the 
principal part of the Laurent expansion of Wfe/(1 —psWk) 
in the vicinity of the point k — —iV{l/s)/li — > 0. Taking 
the next term of this expansion, [ilik+c]~^ (c > 0), which 
is responsible for the behavior of T~^{wk/{1 — PsWk)} 
at hx < 0, one gets J^-^{wk/{l - p^Wk)} ~ e^'^/'Vl^il 
{xli < 0). However, since for small s the condition 
c/V{l/s) >• 1 holds, J^^^{wk/{1 — PsWk)} as a func- 
tion of a; at hx > varies much slower than it does for 
hx < 0. Therefore, in the limit s — > the contribution 
coming from the region with hx < is negligible, and 
the inverse Fourier transform in Eq. (12. 8[) can be safely 
approximated by Eq. p.9p (see also Sec. IIII A| . 

For the same reason as before, in the biased case the 
first term in the right-hand side of Eq. (|2.8p can also be 
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neglected in the limit s — 0. Therefore, using Eqs. (|2.8I) . 
and p.9p . one gets 



6'(Zia;) 



(s — ?> 0), and the Tauberian theorem yields 



-xV(t)/l 



^e{hx) 



(3.10) 



(3.11) 



{t oo). Accordingly, the limiting distribution of the 
scaled particle position Y{t) — X{t)V{t)/li is described 
by the asymmetric one-sided exponential density 

V{y)=\iu, ilipf A^y,A =e-^%). (3.12) 



t-^cx^Vit) \V{t) 



Equations p.Sp and (13.121) are the main results of this 
paper. They show that in the case of the decoupled 
CTRWs with superheavy-tailed distributions of waiting 
times the limiting distribution of the properly scaled par- 
ticle position is described by either the two-sided or the 
one-sided exponential density. Remarkably, these two- 
and one-sided exponential densities correspond to all un- 
biased and all biased CTRWs, respectively. 



A. Testing jump density 



The limiting probability densities (|3.8p and (|3.12p were 
obtained from the approximation of the inverse Fourier 
transform in Eq. (|2.8p by the asymptotic formulas p.4p 
and (|3.9|) . Although the applicability of this approxima- 
tion is well-grounded, it is desirable to check the validity 
of Eqs. (13. Sp and (|3.12l) by using an exact expression for 
that transform. To this end, we consider the two-sided 
exponential jump density 



w{x) = 



e-^'*+, a; > 
e^''-, x < 



(3.13) 



(k-i- > 0). Taking into account that the Fourier transform 
of this density is given by 



Wk = 



[k -f iK,+)[k — IK-) ' 



we obtain 



(3.14) 



(3.15) 



where 



1 



K± = "^2^1^- ' '^+) + 2 - i^+Y +4(1 -P5)k+k_. 

(3.16) 

Since the right-hand sides of Eqs. (|3.14l) and (|3.15l) are 
similar, the inverse Fourier transform of Wfe/(1 — PsWk) 
is of the form of w{x); that is. 



Wk 



K+K- 



l-PsWk) k+ + fi^ le^*^-, a; < 0, 



(3.17) 



and so Eq. (|2.8p in the reference case reads 

Ps{x) 



1 -Ps ^ , (1 -Ps)Ps K+K- 

d[x) -I — 



If 



j 6-^'^+, X>0 

\ e"^^- , X <Q. 
K then Eq. (|XTa yields 



(3.18) 



the corresponding CTRW is unbiased). 



= (i.e.. 
In this case 

— K- = K\/i~ps ^ kV"'''(L/s), <2 = and 
Eq. (|3TT8l) at s -> reduces to Eq. (|375|) . Thus, the hm- 
iting probability density which corresponds to Eq. p.l8p 
with K4. — K_ is given by Eq. (13.81) . as it should be for the 
unbiased case. If k+ 7^ k_ then the biased CTRW with 
li — (k_ — k+)/k^k^ occurs. According to Eq. p.l6p . 
in this case we have 

k± ~ K+K-\ii\e{Tii) + ^v(^je{±h) (3.19) 



as s — 7> 0. Using this result and Eq. p.l8p . it is not 
difficult to find the asymptotic formula for Ps{x) which, 
by applying the Tauberian theorem, gives 



Pix,t) 



K+K^\h\V{t) 



^xV{t)/h^ > 



K+K-\h\'^ + V{t) l^e^i^+'^-h^ l^x<0 

(3.20) 

as t — ?► cxD. While at hx > the asymptotic formu- 
las (|3.20p and p. lip are asymptotically equivalent, at 
hx < they are different. However, as it was mentioned 
earlier, this difference does not affect the long-time distri- 
bution of the scaled particle position ^(i) = X{t)V{t)/li. 
Indeed, taking into account that the exponential term 
^yK+K-ijv{t) at ?/ < tends to zero as t — )■ 00, the limit- 
ing probability density V{y) associated with probability 
density P{x, t) from Eq. is given by Eq. ((XT^ . 



B. Illustrative examples of waiting-time densities 

The long-time behavior of the exceedance probability 
V(t) is the most important characteristic of the CTRWs 
with superheavy-tailed distributions of waiting times. As 
follows from Eqs. p.6p and p. Ill) , it is this behavior that 
is responsible for the long-time behavior of the proba- 
bility density P{x,t). Below we consider three illustra- 
tive examples of the waiting-time probability density p(r) 
which lead to different asymptotic formulas for V{t). 

As a first example, we consider the probability density 

qb'/^ exp[-61n^(77 + T)] 
= T.^. /„ u^^,. — — , (3-21) 



r(l/g,&ln«,/) 



r] + T 



where r{a,x) = dye^^y'^^^ is the upper incomplete 
gamma function [r(a, 0) = r(a)], and the conditions ry > 
1, 6 > and < q < are assumed to hold. The 
conditions 77 > 1, & > and q > are responsible for the 
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positivity and normalization of p(t), and the condition 
q < 1 guarantees that /i(t) is a slowly varying function. 
The last inequality follows from that the limit 



h{^J.T) 



T-i-oo n.(Tj 



= lim exp(— ^gln/iln"^ ^ t) 



(3.22) 



equals 1 for all /i > 0, i.e., the function ft,(r) is slowly 
varying at infinity, only if 5 < 1. According to the defi- 
nition (12.51) , the exceedance probability that corresponds 
to the probability density p. 211) reads 



V{t) = 



r[l/q,b\n\r] + t)] 
r(l/g,&ln«,7) ■ 



(3.23) 



Therefore, using the asymptotic formula [221 r(a,x) 
> 00, or 

51/9-1 



e ^cc" ^ as a; — )• 00, one gets in the long-time limit 



Vit) 



■exp{-b\nH)\n^-H. (3.24) 



Our second example is the probability density 



p{t) 



{r — 1) In*^ ^ r/ 
{i] + t) ln''(?7 + t) ' 



(3.25) 



where r > 1 and i] > 1. The main feature of this density 
is that its right tail is heavier than in the previous case; 
that is, p{t) at t — >■ 00 tends to zero slower than that in 
Eq. p.2ip . The exceedance probability associated with 
the probability density p.25p has the form 



Vit) 



In 77 



ln(77 + t) 



and so 



V{t) 



In 77 
InT 



(3.26) 



(3.27) 



as i — > 00. Comparing the asymptotic formulas (I3.24p 
and (|3.27p . we conclude that the exceedance probability 
in Eq. (|3.26|) decreases slower than the exceedance prob- 
ability in Eq. i^^ . 

Finally, as a third example, we consider the waiting- 
time probability density 



p{t) 



(r- 1) (In In 77)" 



(?7 + t) ln(?7 -f T)[lnln(77 + t)Y 



(3.28) 



(rj > e, r > 1), whose right tail is heavier than for pir) 
from Eq. (I3.25p . In this case the exceedance probability 
is also calculated exactly, yielding 



V{t) 



In In 77 
In ln(77 + t) 



and 



Vit) 



In In 77 
In Int 



r-l 



(3.29) 



(3.30) 



as t — cxD. We note that the above examples are illustra- 
tive and they do not exhaust all possible behaviors of the 
exceedance probability Vit) at long times. Moreover, a 
class of slowly varying functions Vit) contains infinitely 
many functions that, at t — >■ 00, decrease both faster and 
slower than a given function from this class. 



IV. MOMENTS OF THE PARTICLE POSITION 

The moments M„(t) in — 1,2,...) of the particle po- 
sition are defined in the usual way: 



M^it) 



dxx"^Pix,t). 



(4.1) 



Their behavior at long times can be easily determined by 
using the limiting probability densities p.Sp and (|3.12p . 
Specifically, in the case of unbiased CTRWs for the even 
moments we obtain 



M2nit) 



J \ n+1/2 poo I i 7 



2Vit) 



2Vit) 



dyy'"e 



dyy'"P 



2n-y 



2Vit) 



y,t 



(4.2) 



Taking into account that dyy^'^e ^ = r(27t + 1) 
(27i)!, the above result becomes 



^-W-(^)"(277)!. 



(4.3) 



Since in this case Pix, t) is a symmetric function of a;, all 
odd moments equal zero: Af2n-i(0 = 0- 

For the biased CTRWs Eqs. (jH]) and ([332]) give 



- im) 



lioo 



Therefore, since S^2i^oo dyy"'^ ^(^iy) 
asymptotic formula (14. 4p reduces to 



M„(i) 



h 
Vit) 



/i77,!/|/i|, the 



(4.5) 



An important feature of the moments of the particle 
position is that they are slowly varying functions at in- 
finity. This follows from the fact ^20(1 that if some func- 
tion [in our case Vit)] is slowly varying then so does 
any power of this function. Due to this feature, in the 
long-time limit the condition Mnit)/tP -> holds for all 
p > 20]. It shows that the moments Ajf„(i) grow to 
infinity (because Vit) tends to zero as t — > 00) slower 
than any positive power of time. Hence, the CTRW with 
superheavy-tailed distributions of waiting times can be 
viewed as a generic model of superslow diffusion. 
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A. Laws of superslow diffusion 

The character of diffusion is usually determined by the 
law of diffusion; that is, the long-time behavior of the 
mean-square displacement or variance 



a\t)=M2{t)-M^{t). 



(4.6) 



Using the asymptotic formulas (|4.3I) and (|4.5|) . from the 
definition (14.61) we find the diffusion laws 



a' it) 



h/Vit), h^O 



(4.7) 



for both unbiased (h = 0) and biased (h ^ 0) CTRWs 
with superheavy-tailed distributions of waiting times. 
Since in this case the variance cr^(i) is a slowly vary- 
ing function, it grows to infinity slower than any positive 
power of t, i.e., superslow diffusion occurs. As is clear 
from Eq. (14. 7|) . the biased superslow diffusion is faster 
than the unbiased one. We note that the same laws of 
superslow diffusion have also been derived by using an- 
other approach , which does not involve explicitly the 
long-time behavior of the probability density P(x,t). 

Although superslow diffusion has been observed in very 
different systems, all previously known diffusion laws 
are given by a power function of the logarithm of time: 
cr^(t) cx In^i (i/ > 0). The best known example of this 
type of diffusion is the Sinai diffusion [2^ for which v = A. 
Other examples have been observed, e.g., in resistor net- 
works two-dimensional lattices @,|26], CTRWs [11], 
charged polymers [13] , aperiodic environments 2§\ , iter- 
ated ma ps |29j . Langevin dynamics [soj . and fractional 
kinetics [31[. In contrast, the asymptotic formulas (|4.7p 
define a large class of laws of superslow diffusion, which 
includes (j'^{t) cx In" f as a very particular case. This 
class consists of the slowly varying functions that tend to 
infinity as time evolves. 

According to Eq. (14.71) . the diffusion laws are com- 
pletely determined by the asymptotic behavior of the 
exceedance probability V(t). In particular, for the ex- 
amples considered in Sec. IIII Bl the laws of superslow 
diffusion can be obtained by the direct substitution of 
the asymptotic formulas ((XTTI) and ((X^ into 

Eq. (|4.7p . For the first example we obtain 




61/9-1 
/2r(l/g,&ln«77) e 



2 ^2()ln«t 



In'^'-"^ t 

b In" t 



61/9-1 



h = 0. 



(4.8) 



The second example leads to the following law of super- 
slow diffusion: 



In 



2(r-l) 



1 r-1 

L In 77 



t, h^O 



h=0, 



(4.9) 



which for =0 was first derived in Ref. [l3|. In this 
example, that reproduces the diffusion laws with a'^{t) oc 
In"^ t, diffusion occurs slower than in the previous one. 
Finally, the third example yields 




(In In 77)' 



— (lnlnt)2(--i), 1,^0 



h 



(In ln77)'^ 

(4.10) 

i.e., the variance grows slower than that in Eq. (|4.9p . 
We recall that these examples are illustrative because 
infinitely many different laws of superslow diffusion exist. 



1. Specific application 

Now, as an application of our results, we will show that 
particles moving under a constant force /(> 0) through 
a randomly layered medium can exhibit the biased su- 
perslow diffusion. According to [11] , the influence of this 
medium on the particle dynamics can be described by 
the piecewise constant random force g{x) = <?'•"■', where 
X G [nl, {n + l)l), 7T, = 0, 1, . . ., Hs the layer thickness, and 
f;^"-* are independent random variables distributed with 
the probability density u{g) in the interval [—go, go]. In 
the overdamped regime the particle position xt (xq = 0) 
is governed by the equation of motion lyxt — f + g{xt) 
with 1/ being the damping coefficient. At / > go the 
long-time solution of this equation can be approximated 
by the walker position X(t) if the jump density is given 
by w{x) — 6{x—l) and the waiting-time density p{t) is as- 
sociated with the probability density of the random time 
T*^") that a particle s pen ds moving from the point nl to 
the point {n+l)l Since r(") = vl I {f + g^''^) , the 

probability density u{g) which characterizes the medium 
disorder is expressed through p(r) as follows: 



<5) = 



vl 



(4.11) 



It should be noted that the waiting-time density p{t) is 

Tmin, Tmax], whcrC T^in = vl / (/ + 



an) 



defined in the interval 
go) and T,„ax = ^l/U - 

It was shown in Ref. [33| that the biased diffusion of 
particles can be anomalous only \i f — g^. In this case 
the character of anomalous diffusion is determined by the 
behavior of p{t) as t cxd (we recall that Tmax = cxd at 
f = 9o) or, as it follows from Eq. (|4.1ip . by the behavior 
of u{g) as g ^ —go + 0. Thus, taking into account the 
asymptotic formula (|2.9p . one can see that if 



1 



90 + g \9Q+9 



1 



(4.12) 



(g — > —go + 0) then the biased diffusion of particles mov- 
ing under a force f = go a randomly layered medium 
is superslow with (T^(i) ^ l'^/V'^{t) as t — )■ 00 Jsince 



7 



w{x) = S{x — l), we have h = I]- In particular, Eq. (|4.12p 
for p{t) from Eq. ([5:^ yields 

ui9)- . ^r']'f~l\lr ' (4.13) 
[go + 9)\M9o + 9)\ 

and so the law of superslow diffusion in such a medium 
is given by Eq. with = ^ ^ 0. 



V. CONCLUSIONS 

We have studied the long-time behavior of the de- 
coupled CTRW with a superheavy-tailed distribution of 
waiting times. We assume that the distribution func- 
tion of waiting times belongs to a class of slowly varying 
functions and the jump distribution has a finite second 
moment. At these conditions, the limiting probability 
density of the properly scaled particle position has been 
derived for both unbiased and biased CTRWs. In the for- 
mer case the limiting density is given by the symmetric 



two-sided exponential density, while in the latter case it 
takes the form of the asymmetric one-sided exponential 
density. We have used these limiting densities to calcu- 
late the moments of the particle position at long times. 
It has been shown that all the moments are expressed 
through the exceedance probability and grow slower than 
any positive power of time. 

Due to the last property, the decoupled CTRW with 
a superheavy-tailed distribution of waiting times can be 
viewed as a generic model of superslow diffusion. We 
have derived the most general form of the laws of su- 
perslow diffusion and have shown that, in the cases of 
unbiased and biased diffusion, the variance of the parti- 
cle position is inversely proportional to the first and sec- 
ond powers of the exceedance probability, respectively. 
The laws of superslow diffusion that correspond to the 
illustrative examples of superheavy-tailed distributions 
of waiting times have also been determined. Finally, we 
have applied the obtained results to show that the biased 
diffusion of particles moving under a constant force in a 
randomly layered medium can be superslow. 
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